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PREFACE 

The research described here was performed for the National 
Aeronautics and Space Administration, and deals with the APOLLO 
Mission Reliability Assessment Study. In this Memorandum, the 
author uses Bayesian analysis to specify parameters of a prior 
distribution for two cases: (1) reliability of a unit that 

either performs satisfactorily throughout a mission or does 
not, and (2) failure rate of a unit that fails according to 
the exponential distribution. Prediction of demand for spares 
is considered in each case. The cases can be read independently. 

An estimate of reliability is the posterior mean. Alterna- 
tively, the posterior distribution can be used to obtain a (sub- 
jective) confidence interval for reliability. The posterior 
distribution is also useful in a decision- theoretic approach 
to resource allocation for maximal system reliability; such a 
study is planned as a sequel to the present work. 

This Memorandum should be of interest to those working on 
reliability estimation; allocation of investment among system 
components to achieve maximum system reliability; and stockage 
applications . 
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SUMMARY 


This Memorandum specifies the parameters of a prior distribution 
for two cases: the reliability of a unit that either performs satis- 

factorily throughout a mission or does not; and the failure rate of a 
unit that fails according to the exponential distribution, 

Bayesian analysis is an obvious approach in estimating reliability 
parameters from mixed data sources such as: (1) test results; (2) 

information on analogous components; and (3) engineering estimates . 

The prior distribution, of necessity subjective, is (ideally) based 
on (2) and (3) alone. Test results are then merged with the prior 
via Bayes' rule to obtain a posterior distribution. Roughly, the 
spread of the prior distribution is inversely proportional to the 
degree of prior belief, and determines how heavily it will be weighted 
when combining it with test data. 

A topic that most writers on Bayesian analysis avoid is how to 
specify the parameters of the prior distribution based on (2) and (3). 
We give a method for specifying these parameters that requires only 
information corresponding (i) to the most likely value of reliability 
and (ii) to the subjective odds that the error in this estimate is 
less than a given percent. Wp have computed tables of parameters of 
the prior distribution corresponding to these subjective inputs. 

These appear in the Appendix, 

As an application of Bayesian analysis, we consider prediction 
of demand for spares in both the GO NO-GO and constant failure rate 


cases. 
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1, SPECIFYING THE PARAMETERS OF A PRIOR DISTRIBUTION FOR 
RELIABILITY OF A GO NO-GO UNIT 


Suppose we have a unit that works with probability p but that 
the precise value of p is unknown. If we were totally Ignorant about 
the value of p, then our prior belief would be reflected by a uniform 
distribution over [0,1]. However, Intuition tells us that total 
ignorance is an anomaly; that is, our prior distribution is really 
not flat. A smooth, unimodal prior distribution having support (0,1) 
with peak over what we believe to be the most likely value of p seems 
appropriate. In addition, the beta distribution is a natural con- 
jugate [5] prior distribution; I.e., the posterior distribution is 
again a beta distribution (with parameters transformed according to 
Bayes' rule). The beta density with positive parameters (a,b) is 
given by 


( c p a x (l-p) b 1 , 0 < p < 1 

(1.1) P(p|a,b) = ] 

( 0, elsewhere 

with c as a normalizing factor. The mean and variance are, respectively: 
(1*2) ^ = a/ (a+b) , 

(1*3) = ab/[(a + b)^(a + b+l)], 

and, for a, b > 1, there is a unique mode at 
(1.4) 0 = (a- 1) / (a+b-2) . 

After observing test data, say a sample with m successes and 
n failures, the posterior density is p(p|a+m, b+n) from the Bayes' 
rule relation: posterior density * prior density x likelihood function 

x a normalizing factor independent of p. As more test data are observed, 
the posterior distribution is updated. (The updating procedure Is valid 
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only if all units are stochastically identical- If, for example, 
design changes are made, as a result of failure mode analysis, a 
new prior distribution should be constructed from scratch) - 

It remains to specify the parameters (a,b)- The procedure we 
give is heuristic and, while not the simplest mathematically, uses 
information that corresponds to subjective notions. For example, 
one is less likely to have a feel for the variance of the prior 
distribution than for the error in his estimate of the most likely 
value of p. Of course, if we were interested in psychological con- 
sistency, we could ask for an estimate of the variance as well -- 
but we shall ignore such considerations here. If the designer/ 
engineer being asked these questions has seen any test data, it is 
probably impossible for him to ignore them. Therefore, in this 
case, it is suggested that the prior distribution be based on all 
information the designer knows. On the other hand, if the designer 
has not seen any test data, this is all to the good; test results 
are then accounted for in the posterior distribution. Whenever 
possible, the parameters of the prior distribution should be speci- 
fied before any tests are performed. 

Suppose that our subjective assessment of the most likely value"* 
of p is p; then we set 

(1*5) (a- 1) /(a+b -2) = p . 


For example, if £ and £ were subjective estimates of the mean 
and variance, respectively, of the prior distribution, then solving 
the equations (1.2) and (1.3) yields 

a = o £ 2 (i-£) - (1 
b = a" 2 jl(l-|i) 2 - (l-(i). 

icic 

The analysis of the case where one estimates the mean rather 
than the mode is analogous. We give no details for the former case, 
except that numerical results for both cases are given in Tables 1 
and 2 in the Appendix. 
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Next, we ask what odds we would give that the true value of p lies 
in (p-kp, p+kp) , where 0 < k < 1. For example, if k * .1, then 
we ask what the chance is that the error in our estimate is less 
than 10 percent. Suppose that the subjective odds are x to y; 
then, setting v = x/(x+y) , we have 

J$+kp 

(1.6) I p(p|a,b) dp = v. 

p-kp 


Thus, to find a and £, it suffices to specify p, k, and v. If 
the views of several people are solicited, it is suggested that the 
decision-maker take weighted averages, the weights {cy^} depending 
on the technical backgrounds and personalities of the people in- 
volved. Some may be conservative in their estimates, while others 
are optimistic. It is suggested that in asking the questions the 
decision-maker fix the value of v. If person i responds 
then p - la^Pj/Ecfj and k = 

Equations (1.5) and (1.6) can be resolved by using the tables 
of the incomplete beta function [3], but to expedite matters we 
have provided a table of (a,b) in the Appendix corresponding to 
selected values of (u,v,k), where u ■ f>. 

Defining 


•min(l ,u(l+k)) 


(1.7) 0(t;u,v,k) = v 


iu( 1- 


(1-k) 


u(t-l) 

1-U . .t-1 

P (1-p) dp 


1 u( t -U 

f p l '“ (i-P ) 1 ' 1 

* / 0 


dp 


and 


( 1 . 8 ) 


0[g(u,v,k) ;u,v,k] *= 0 , 
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it follows from (1.5) and (1.6) that 
(1*9) b * g(p,v,k) , 

(1.10) a * [p(b-2) + l]/(l-p) . 

A uniform prior is suggested, if there is not enough prior information 
to quantify sensibly; however, it is felt that introspection will 
generally reveal the contrary. 


In base stockage application [2], appropriate levels of spares 

inventory must be determined. To provision spares properly, an 

estimate of the demand distribution is required. For this, the 

Poisson approximation may be useful. Assuming that p is near one 

and the sample size n (say, aircraft landings or space vehicle 

* 

launchings) is large, the probability of k failures, corresponding 
to demands for spares of a given type, is closely approximated by 

(1.11) f(kjp,n) = [n(l-p)] k e' n(1_p) /k:. 

Removing the conditioning on p, the demand distribution is 

(1.12) g(k|a,b,n) = I f(k|p,n) p(p|a,b) dp. 

J 0 

An approximation to g(k|a,b,n) is obtained by using the mean of the 
prior distribution [a/(a+b)] in place of p in (1.11). We do not 
know how good this approximation is. 


We assume that the failure distribution over successive missions 
is geometric (i.e., no memory). The part in question is assumed stressed 
(used) exactly once per mission -- or, with obvious modifications, twice 
per mission. If it is stressed continuously, the results of Sec. 2 can 
be applied; of course, if all missions have the same length, we get a 
reduction back to the case considered here. 
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If the distribution of n, say 0(n) , is known, then the distri- 
bution of the number of failures is 

h (k | a , b) = ^^g(k|a,b,n) 0(n) . 
n 

To the author this indirect route to demand prediction seems preferable 
to a direct attack because the former is more physically motivated. 

A device sometimes used is to inflate the estimate of demand 
deliberately in order to cause a larger provisioning of buffer stocks, 
with the object of reducing the incidence of stockouts due to demand 
fluctuation. The author feels that the approach outlined in the next 
paragraph is more rational. 

With an unbiased estimate of demand, the proper Inventory level 
can be determined by a decision- theoretic approach. Let L(k,s) be the 

’ A * ' A * 

loss when k units are demanded and s units are stocked. The minimal 
expected loss L° is 


(1-14) L° = min L(k,s) h(k|a,b) . 

k=0 


Minimizing L(k,s), where k is the estimate of mean demand, may be 
grossly incorrect. 


A- 

Predicting n via a Bayesian approach -- perhaps in conjunction 
with spectral analysis of time series -- may be appropriate. This 
involves simply one more conditioning-unconditioning in (1.13). Since 
the values of n over successive time periods may be autocorrelated , 
spectral analysis may be useful in finding a suitable parametric form 
for the distribution of n. For a treatment of spectral analysis, see 
Yaglom [6]. 

** 

For example, 


L(k,s) = c^ s + c^ max(0,s-k) + c^ max(0,k-s) , 

where c^ is the unit purchase cost, c 2 is the unit holding cost, and 
c^ is the unit stockout cost. 
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2. SPECIFYING THE PARAMETERS OF A PRIOR DISTRIBUTION FOR 

FAILURE RATE 


Suppose that a unit has constant failure rate 0, fixed but 
unknown. We assume a (natural conjugate [5]) prior distribution 
with density of the form 

(2-1) g( 0 |a,b) = a b 0 b "" 1 e" a9 /r(b) , 

where the parameters (a,b) are to be specified. Its mean and variance 
are: 

(2.2) = b/a , 

(2.3) o ^ m b/a^, 

respectively. There is a unique mode at (b-l)/a, b ^ 1, but it is 
felt that the mean time to failure is more amenable to subjective 
assessment in this case. 

If the subjective estimate of the mean time to failure is v, 
then using (2.2) we set 

(2.4) b/a = l/ v . 

Based on subjective odds, let v be the chance that the failure rate 
exceeds k/\S. This yields 


(2.5) 


| g(eja,b) de - V, 

k/o 
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Equations (2.4) and (2.5) could be resolved by using tables of the 

ic 

incomplete gamma function [4], but this would be a tedious job. 

To save time, for selected values of k and v, Table 3 of the Appendix 

A 

provides the corresponding h, where 


(2.6) 

a - hv, 

(2.7) 

A A 

b = h. 

Defining 


** 

6 (h ;k ,v) = v - I g(e|h,h) d0 , 

(2.8) 



(2.9) 

6[p(k,v) ;k,v] = 0 , 


we see that 


* 2 
If we had used a subjective estimate, say a , of the variance 

of the prior distribution instead of (2.5), then using (2.3) and (2.4) 

we would have the explicit expressions 

a = l/va 2 
b = l/(0a) 2 . 

2 

However, it is felt that intuition for a would be poor. 

’k'k 

Note that 



h,h) de = r(h,hk)/r(h,0) 




where 


r(a,x) 



du. 


A standard subroutine for computing f(a,x) is available. 
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(2.10) h = p(k,v). 

In specifying the parameters of the prior distribution, we refer 

to the suggestions given in Sec, 1 for handling data already on hand. 

Failure data (except that used In forming the prior distribution) are 

incorporated in the posterior distribution by Bayes’ rule. Having 

observed failures at ages t , . .., t , and a nonfailed group with 

1 R * 

ages t , . t , the posterior density is 


g(e| a + ) V ,b+k). 


i=l 


This gives us our current prior distribution, which is updated as more 
observations are recorded. If, for example, the unit corresponding to 
t , j > k, fails at t ', updating yields 

m 

i = l 

We now give an application to demand prediction. Suppose each of 
n units operates continuously until failure, at which time it is 
replaced instantaneously (for practical purposes) by a unit as good 
as new. These failures generate the demands for spares and/or repair. 
If each unit has constant failure rate 0, the probability of k demands 
in time T is p(k|n0T), where 


( 2 . 11 ) 


P(k|\) = X R e" Xk /k!; 


•flf — Qy 

If the failure distribution were 1-e 

O' 

shape parameter a) > replace t^ by t^, i ■ 1, 


a 


(Weibull with known 
• , m. 
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that is, demand is Poisson with rate n 0 . Removing the conditioning 
on 0 , the probability of k demands is 


( 2 . 12 ) 


f (k|a,b,n,T) 



p(k|n 0 T) g( 0 |a,b) d 0 . 


If n and/or T is a random variable, we can further uncondition in a 
similar manner. 


REMARKS. 


y- 


j 


and is 


If, in fact, the life distribution of the jth unit has mean 
nonlattice but not necessarily exponential, then [ 1 ], with 


e = 


n 

< i/n >L 


j=i 



and n-»oo , 


the stationary demand distribution becomes p(k]n 0 T). If planned 
replacement takes place at age t, the same result holds if we replace 
by the mean of the distribution truncated at 7 . (In a more sophisti- 
cated replacement policy, the planned replacement age should ideally 
depend on the current inventory level.) If replacement can take a 
significant amount of time (due, for example, to stockouts or non- 
negligible repair times), then the replacement time distribution should 
be convolved with the failure distribution, and the mean of the result- 
ing distribution used in place of j* . 


For aircraft spares provisioning, a somewhat different model of 

the demand process may be appropriate* Suppose that a part, used only 

when the aircraft is flying, has constant failure rate 0 during flight 

and failure rate 0 on the ground. Let flights to the base originate 

from points fl, ml, with respective flying times fw n , w 1 . 

^ 1 


See Sec. 1 for the case where the unit is not stressed continu- 
ously during flight. 
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Durlng a period of length T, let be the number of flights to the 

base from point i. The probability of k demands during time T is 


(2.13) D(k j T) = 


1. \ o ' ' 


(k l‘"’ ,k m )eS k 


where 


(2.14) 


0w. 


v . * 1-e 

i 


m 


(2.15) 


\ ■ { <k i ■ k } 

i=l 


If all the v.'s are near 0 and all the n.'s are large, then we have 
l i 

the Poisson approximation 


(2.16) 


D(k|T) 


« P 




When the n^'s and 0 are unknown, we condition and then uncondition 
in the usual way. If 0 has a prior distribution g(0|a,b) given by 
(2.1), then 


( k ii>iV a ' b ) ■ 4e "if 1 '(^i) ] ) 


(2.17) 
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AP PEND IX 


In this Appendix we give tables of parameters of prior distri- 
butions corresponding to subjective assessments of reliability, as 
described in Secs. 1 and 2. The entries in the tables were computed 
using a program written by Mrs. Sarah Higgins, with the assistance of 
Robert Mobley and the author. Several test cases were computed by 
hand (using tables) for each program, with agreement to more than 
four significant figures throughout. The programs are believed to be 
completely debugged and are listed here for the convenience of those 
who may want values that are not tabulated. Tables 1 and 2 refer to 
Sec. 1 (beta prior). Table 3 refers to Sec. 2 (gamma prior). Aster- 
isks in the tables indicate that the rootfinder did not locate a root 
within the allotted time. 
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Table 1 


PARAMETERS OF BETA PRIOR DISTRIBUTION FOR SELECTED VALUES 

OF u, v, AND k 


A Priori 
Estimates 

Mode 

Mean 

u 

■91 

k 

a 

b 

a 

b 

.9 

.5 

. 1 

5.510 

1.501 

15.311 

1.701 



.075 

8.995 

1.888 

18.856 

2.095 



.050 

19.067 

3.007 

29.012 

3.224 



.025 

73.625 

9.069 

83.657 

9.295 


.67 

.1 

11.233 

2.137 

21.405 

2.378 



.075 

18.352 

2.928 

28.576 

3.175 



.050 

39.311 

5.257 

49.576 

5.508 



.025 

* 

* 

163.424 

18. 158 


.75 

. 1 

15.862 

2.651 

26.222 

2.914 



.075 

25.727 

3.748 

36.173 

4.019 



.050 

54.915 

6.991 

65.395 

7.266 



.025 

* 

* 

if 

if 


.95 

.1 

50. 195 

6.466 

61.012 

6.779 



.075 

80. 184 

9.798 

91.500 

10. 167 



.050 

165.106 

19.234 

176.772 

19.641 



.025 

* 

* 

if 

if 


.99 

. 1 

90.232 

10.915 

101.161 

11.240 



.075 

145.467 

17.052 

157.047 

17.450 



.050 

* 

* 

if 

if 



.025 

* 

if 

if 

if 

.95 

.5 

. 1 

5.910 

1.258 

25.348 

1.334 



.075 

7.954 

1.366 

27.555 

1.450 



.050 

12.165 

1.588 

32.062 

1.687 



.025 

39.270 

3.014 

59.320 

3.122 


.67 

.1 

10.828 

1.517 

30.424 

1.601 



.075 

15.145 

1.744 

34.972 

1.841 



.050 

24.645 

2.244 

44.926 

2.365 



.025 

80.430 

5.181 

100.846 

5.308 


.75 

.1 

14.503 

1.711 

34.177 

1.799 



.075 

20.667 

2.035 

40.602 

2.137 



.050 

34.680 

2.773 

55.141 

2.902 



.025 

112.251 

6.855 

132.911 

6.995 


.95 

.1 

39.402 

3.021 

59.307 

3.121 



.075 

59.151 

4.061 

79.382 

4.178 



.050 

107.748 

6.618 

128.632 

6.770 



.025 

339.228 

18.801 

361.167 

19.009 


.99 

. 1 

67.365 

4.493 

87.362 

4.598 



.075 

102.944 

6.366 

123.272 

6.488 



.050 

192.310 

11.069 

213.312 

11.227 



.025 

* 

* 

* 

* 

.975 

.5 

. 1 

6.239 

1.134 

45.458 

1. 166 



.075 

8.411 

1.190 

47.710 

1.223 



.050 

12.734 

1.301 

52.197 

1.338 



.025 

25.695 

1.633 

65.637 

1.683 


.67 

.1 

10.689 

1.248 

49.979 

1.282 



.075 

14.711 

1.352 

54.111 

1.387 



.050 

23.166 

1.568 

62.790 

1.610 



.025 

51.772 

2.302 

92.096 

2.361 
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Table 1 -- continued 


A Priori 
Estimates 

Mode 

Mean 

u 

IB 


a 

b 

a 

b 

.975 

.75 

.1 

13.847 

1.329 

53.174 

1.363 



.075 

19.258 

1.468 

58.709 

1.505 



.050 

30.927 

1.767 

70.631 

1.811 



.025 

72.607 

2.836 

113.105 

2.900 


.95 

. 1 

33.943 

1.845 

73.393 

1.882 



.075 

48.848 

2.227 

88.461 

2.268 



.050 

83.263 

3.109 

123.202 

3.159 



.025 

223.049 

6.694 

263.966 

6.768 


.99 

. 1 

55.718 

2.403 

147.780 

2.442 



.075 

81.364 

3.061 

167.403 

3.104 



.050 

141.906 

4.613 

181.930 

4.665 



.025 

396.793 

11.149 

437.834 

11.227 

.99 

.5 

.1 

6.441 

1.055 

105.529 

1.066 



.075 

8.695 

1.078 

107.815 

1.089 



.050 

13.190 

1.123 

112.374 

1.135 



.025 

26.601 

1.259 

125.981 

1.273 


.67 

.1 

10.592 

1.097 

109.706 

1.108 



.075 

14.423 

1.136 

113.580 

1. 147 



.050 

* 

* 

121.496 

1.227 



.025 

47.093 

1.466 

146.599 

1.481 


.75 

.1 

13.437 

1.126 

112.565 

1.137 



.075 

18.381 

1.176 

117.556 

1.187 



.050 

28.610 

1.279 

127.882 

1.292 



.025 

62.041 

1.617 

161.611 

1.632 


.95 

. 1 

30.642 

1.299 

129.821 

1.311 



.075 

42.607 

1.420 

141.851 

1.433 



.050 

68.396 

1.681 

167.770 

1.695 



.025 

160.581 

2.612 

260.346 

2.630 


.99 

.1 

48.573 

1.480 

147.780 

1.493 



.075 

68.123 

1.678 

167.403 

1.691 



.050 

111.020 

2.111 

210.442 

2.126 



.025 

269.774 

3.715 

369.615 

3.733 

.999 

.5 

.1 

6.565 

1.006 

1005.572 

1.007 



.075 

8.871 

1.008 

1007.881 

1.009 



.050 

* 

* 

1012.496 

1.014 



.025 

27.293 

1.026 

* 

* 


.67 

. 1 

10.530 

1.010 

* 

* 



.075 

14.241 

1.013 

1013.255 

1.014 



.050 

* 

* 

1020.678 

1.022 



.025 

44.127 

1.043 

1043.174 

1.044 


.75 

.1 

13.186 

1.012 

* 

* 



.075 

17.842 

1.017 

1016.857 

1.018 



.050 

* 

* 

* 

* 



.025 

55.493 

1.055 

1054.545 

1.056 
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Table 2 


PARAMETERS OF BETA PRIOR DISTRIBUTION, VARYING u 
(v = .5, k = .2) 


A Priori 
Estimate 

Mode 

Mean 

u 

a 

b 

a 

b 

.990 

3.049 

1.021 

102.089 


.991 

3.055 

1.019 

113.201 

1.028 

.992 

3.060 

1.017 

127.091 

1.025 

.993 

3.066 

1.015 

144.950 


.994 

3.072 

1.013 

168.762 

1.019 

.995 

3.078 

1.010 

202.097 


.996 

3.083 

1.008 

252.099 

1.012 

.997 

3.089 

1.006 

335.434 

1.009 

.998 

3.095 

1.004 

502.103 

1.006 

.999 

3.100 

1.002 

1001.997 



Table 3 


PARAMETERS OF GAMMA PRIOR DISTRIBUTION FOR 
SELECTED VALUES OF k AND v 


k 

m 

A 

h 

. 1 

.5 

* 


.67 

.3683137 


.75 

.4964839 


.90 

.9725192 

.5 

.5 

.5602821 


.67 

1.4046943 


.75 

2.1610065 


.90 

5.3209309 

.9 

.5 

3.2660242 


.67 

* 


.75 

* 


.90 

* 
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? 
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4 
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To 
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?s 
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"70 F »R AT {F10.4 f ?<sX,Fs.'4) f 2<*>X,Kl0.4) t ->* f hlo.»<) pp 

VOO F l iRr-i A T p^ 

uul FM'Jh'jT (SFS.3) 30 

VO? FOR-iAT (4Ps.3) 31 

03 FOR -AT {Fs.l) yj 

R-AO ( 3, HSU ) I I 3 3 

RPAI) (H,RS)) JJ -V, 

KhAI' (3,HS0) << 33 

RPAO (S,VOO) ( 0{ I } , 1 = 1 t I I ) ^ 

RF AO (R,VOl) J = KJJ) * f 

RFAI) |S f Vo?) ( X<(<) * <«1 t <<) ^h 

RFAO (s,V03) TX 3 * 

-'=1 41 

P =? ) 41 

IF = l) 4 ? 

otj BOO K* l ,KK 43 

on svv 1=1,11 44 

M.I S9H 43 

f f ST*n ( I )*< 1 . + XU <) ) 43 

IRPhR*AoI l (1 . , f PSi ) t, ( 

?3 XUl^Ksl!lI)»fU-X<|<n 4 H 

Call Ok T f iMfTXf lH t IF) 4w 

*■ (H( I )*( TX-?. )♦!.)/( I .-.III) ) 3 1 

-Mils (f5 f M7o) H(I]|V(J) 9 X<(K),A v TX v PHl Si 
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Call pxit 
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Hftl'TC G*TF i. 1ST 

C G€NE-*Ai_ KMlTr Ini)*:* V-UTTHN HV wEK-'Hk 

SUMKHUTlMe G«T iMfCfl^tlr) 

OIM6NSIMN C < 50 > 

oo too l* l ■» n 

JK*0 

IF (C(U )45, 46,45 
45 *T«.V*CIU 

ASSIGN 1 TO mn 

GO r»lT HO 
1 XO*F PKT 

*T«1.1*C(U 

ASSIGN 2 TO nn 
GO TO 80 

2 X 1 * F P K T 

>‘T«C(U " 

ASSIGN 3 TO NM 
Gu TO 8o 

3 X2* F PKT 
GO TO SO 

46 RT*- I .0 

assign A TO nn 
GO TO 80 

4 XU*FPkT 
KT-1.0 

ASSIGN 5 TO MW 

GO TO H0_ 

s ”x i »fpkt 

KT*0*O 

ASSIGN 6 TO MM 
GO TO 80 
' 6 X2*FPk T 

_50 H«- 1 ,0 

n*» - # 5 

49 ni)*l # u+n 

8 1 * < xo*n*o)-< x i*i)o*oo) > ( x 2 *(oo«h>) ) 

OF>'*R I*H I- I 4 # 0*XZ*l>*»)0)*< XO*l)-{ Xl*I.M)) + X2 ) 

IF (OEM) 36,36,51 

_ 36 0_FN*Oj, 0 

SI n H n* $ UKT ( 0Fn ) 

53 OM*bI+OEN 
0**8 1-OHM 

IF ( A8SI0N)-ft8S(0w) 157,57,56 

54 OFm*dn 

GO TO 58 

6 7 OHM *04 

*58 IF (0FN)55, 54,55 

54 

55 0I«I-2*0*X2*D0)/06n 
H*0 I *H 

8T«KT + H 

\f < AHS(H/KTT-U0H-4)y5,*75t60 
60 ASSIGN 7 TO nn 

GO TO 80 

7 IF ( A6$ { FP8T )-A3S { X2M0.0 ) ) 62 , 41 , 61 

fei 01 *01 * . 5 

H»H *. 5 

KT-KT-h" 

GO TO 80 
62 X0*X1 


“0)080010 41 

F K AM< - t *C T 03 FR 20, 1958 6? 

■•<0080020 6 3 

i()08003 ii 64 

N0OHO04" 6S 

MOOHOOSn 66 

N0(»8006n 67 

■•'0060071* 6 8 

•^OOrtOOnif 6 v 

*>U0h009o 7 V ) 

'0060100 7 \ 

-00 601 in 1? 

’-'OOhO 120 73 

•'00 *0 1 3 1 74 

0)080140 76 

•’00 801 SO 7 6 

**00 8 0160 7 7 

-'00 801 7* » 1 h 

■<00801 ho 7s 

-•0080 Ivn 80 

■ J 0080?0() Hi 

900*0210 8? 

N00H022O *3 

•*»00H023o 84 

-0)080240 ^S 

•<00«025o 86 

’•'006026c- 6 7 

,,, OQ8027(> 8H 

NQ0HO2HO 

900 HO 2 Vl) VO 

*0080300 vl 

’•'0080310 V7 

NOf)M032f) V3 

’-'0080330 94 

'■'00 8 0340 V8 

■*'U08035o 9 6 

-'0080341? 97 

■■00603 70 Vi 

N00H03KO VV 

■*'0080 3 91! 100 

•10080400 101 

w (JO 6 04 1 0 102 

■■'0080420 103 

*‘■00804 3 O \04 

NU080440 1 OS 

'•i 00 HQ 480 10 4 

WOO 80460 107 

90080470 108 

W(jn«n480 109 

900 8049 0 HO 

60060500 1 1 1 

900H0510 112 

90060520 113 

■900H0 5 3 0 114 

•80080540 US 

■•<0080550 1)4 

MOO 805 60 J17 

90080S70 118 

J 0 0 8 0 b 6 f ; llv 

'■'0080891' 120 
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X 1* X 2 


W0080600 

J?1 

*■ * 

X2-FPKT 


W00B061 0 

1 ? 2 


0«DI 


W0080620 

123 


GO TH 49 


90080630 

f2*~ 

75 

CALL AUX (RT,FRT) 


W0080640 

125 

76 

C(U*KT 


W0080650 

126 

100 

continue 


WOO 80660 

127 


IN»JK 


W00H067O 

128 

33 

RETURN! 


W0080680 

.. .. . 12.9 

HO 

JK«JK+1 


WOO H 0690 

130 


IF (100-J<)75,75,86 


W0080700 

131 

86 

CALL AUX (RT,FRT) 


W00H0710 

132 


FPRT-FRT 


WQOH0720 

133 


IF (L-l)81, 89*81 


W0080730 

134 

81 

nn 8? I ■ 2 • L 


W0080740 

135 


TEM.RT-CI I-I ) 


W008075O 

136 


IF ( A85(TEM}-1 ,0E-20) 85, 82,82 


W00H076U 

137 

82 

FPRT*EPRT/TE™ 


W 0 0 8 0 7 7 0 

138 

89 

IF { AR5(FkT )-l •06-20 >90,91 ,9 1 


WOO BO 780 

139 

90 

IF (A8S(FPRT ) - 1 .Ot-2 0)76,91 ,91 


WOOK(>79i' 

140 

91 

IFIIF) 33t«4 ? 33 


'J 00 80 8 00 

141 

84 

GO TO nn, 11,2,3,4,5,6,7) 


W008081 0 

14? 

85 

RT * RT+ • 00 1 


WOO 80820 

163 

88 

GO Tfl RO 


W008063O 

144 


ENO 


WOO 80 840 

146 





146 





147 





148 





149 

4>I8FTC AUX 



150 

C 

auxiliary program calleo hy koutfi^oek 

SOrtROUT I nE 


151 

C 

this SimnoTHF RETURNS CONTROL T « * THE 

ROOTrlNOFK ROUT I M 

E 

15? 


SIIHKOUTIwe AUXUr^KTJ . . _ ... 



...1*3. . _ 


EXTERNAL FmC 



16 4 


COMmOM/ I NPHT/IH9) , V( 5) ,X<(4) , UPPER, XLDWER t tX,I, 

PHI , J, AX 


155 

598 

FORMAT (4H Wl*,lPE15.7,4H W?*,1PE15.7) 



156 

599 

FORMAT (6H XnOh* ,1 PE 15*71 



157 

600 

FORMAT (5H OEM- , l PE 15. 7 > 



158 

Wi 

-Fjj.^wAj_L4h 



JL12 

60? 

FORMAT ( 4H Z2»,IPE15.7> 



1 60 

HO 3 

FORMAT (5H PH I * , l Pfc 1 5# 7 ) 



161 

806 

FORMAT ( 4H0RT*, 1PE15.7) 



16? 


AX-KT 



163 


OIF* .0001 



164 


WRITE (6,804) RT 



166 


IF (RT-l.) 25,28,35 



164 " 

25 

PHI*V( J )- (OPPFR-XLOVER)*! .-RT 



167 

C 

ARTIFICIAL VALUE iNjFNOtO T» » ORIVE ROUT F 1 AM'gR TOwARO A ROOT 


168 


GO JO 46 



1 69 

28 

PHI «V( J )- (UPPER- XLOwER ) 



170 


.m TO _ .. 



171 

35 

CALL K I N T 2 ( 0 ( I ) , UPPER, 0 IE, FW'’ , R 1 , IN*)) 



17? 


.JIALL 8JL*!X2_ J. x L Q w E r , U ( I ) , OlF,F.v>c t w2 , 1 »m > ) 



173 


WRITE (6,598) Rl,W? 



174 


XNi ih» wl-f w? 



176 


HR J T E (6,599) XNUn 



176 


TOI *0,000 1 4XN0M „ 



177 


X-0.0025 



178 

5 

VALUE«X**(U( I )*(RT-l. )/( l.-UU ) ) ) * I1,-XJ**(RT 

-1. ) 


179 


IF ( VALUE. GE.TOL. ANO* X . FO .0 .002 5 ) GO TO 50 



1 HI) 



onHtnnnonnfvfton ojo n o o 
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IF ( VALUP.LT. TOL*AnO.X. go. 0,002 5) GO TO 1> l \ u \ 

IF ( VALUF.LT .TOC ) Go T»J 39 
X»X LflwiE K-0» 1 

27 I F I X . L E • 6 ♦ 0 ) ’ GO T 0 *>0 lW4 

29 VALUE»X**(IJt I ) *(RT-l. )/ ( l.-U( I) ) ) * (1 .-X ) ** ( KT- 1 . ) 

IF ( VALOF.LT .TO'.) GO TO AO 

X-X-O.l l ^ 7 

r,n to 27 lKH 

32 CALL K I nT 2 I YUWF K» Xjj1 W6H»Olh»FNCf Zi f I -M n ) 

30 w R I | k (*,601) Zl " ' " lv0 

IF (UPPFR.FU. I .0 ) GO TO 40 1^1 

CAt.L RINT? (UPPER, WO.o1E # EnC,72, IND) 1^2 

WRITF (6,602) Z2 

41 0F^*XMUM+ Z l+Z 2 1** 

WRITE (4.600) OEo 

IF ( AbS ( XMU»*» ) *Gr*l.ElO*UEN) GO To 3 ***» 

GO TO 43 * v7 

40 Z2-0.0 * s "' 

GU TO 41 I 4 * 1 * 

43 PH I * V ( J ) — ( X'vU»*i/0«flM ) -00 

GO TO 4 6 " 20}. 

3 PHl«-.001*Ml 202 

46 WRITE (6,hV3) PHI 2^3 

GO TO HO ? (]U 

39 Zi«0.0 

GO If) 30 ?n<7 

3 ) Y LOwHK»0.O. ?r : 7 

GO TO 32” ” ?n>4 

5! X*XL')wPR 2'^ 

<;u to ‘3 21.0 

60 YLOPFK*X 2 1 l 

GU TO 32 ?]? 

HO FRT«P H I 213 

RETURN ' ~ ? } /f 

F iv h 213 

214 

217 

21 9 

SIRFTC R]nT2 LIS1 1037001^ 2?<1 

I n f P GK A T 1 1 mm SUHRI lUT I^'F WRITTEN hY RORErT L. •iO-sL-Y 221. 

iiA T F f}F WR If E-Up - 2-?-6*> 9oS7o'02o 2?2 

n A T P OF SOURCE DECK - 2-2-63 J 0370u3u 223 

SUHKUUT I <MF R I mT 2 ( A ,3 , PHC ,F , I^D) <u67ooao 22S 

•0)37 00 So 2 2*> 

™ *0 37oo 6(i ?2S 

A ■ owe L I i‘i I T OF THE H TEGRAT ION. W037O07O 22/ 

h s other l I ■** I f i if i he Integration. ousyoohc 2?h 

F » ERROR ROUND ( nun-OIREnS IU-mA>. ) . J0S70U9M ?2v 

" FnC * FUNCTION SUKPROGkM. ,J J‘WOl(>u 230 

the function s ( a tf "E *st «e - function f^cixj-'i . *>70110 231 

" ~JjhERE X IS The I Nt)E PFnDF 1 f V 4 K I A h L E • i*io 370) ?o 232 

E « THE VALUE OF THE INTEGRAL IS R- TURNED HfcK E • *103701 3o 233 

I NO * AN I <n»'J ICA TOR RHICH IS RETUROFD. wo37014o 234 

ZERO INDICATES THE IvTEGR*L » * I » » *'OT CunVhPGE '•'O'sYoibU 23-> 

USING ?**10 P'OFkVa'.S# •*•(>6701.60 234 

,viu,m>zeko indicates the integral converged wimi*' <jo ; «/oi7u 237 

♦ • -* PR <1 IK WOUND* -'(n/'llHi. ?Hi 

wo37nJso 23 v 

DOUhLE PRECISION n 31) ,AA»3H,SIGnA f hNC,Ht) f P, E^I.V,R,FF,Au$, A4<,DA 240 
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1HS 

(PM I N» 1 0000 000# 

A A- A 

BH* H ' " " 

HO-HA-AA 

T * • 5 * ( F M C ( A A )+FnC (H3) ) *H0 

P-1. 

DM 100 1*1,1? 

__ S 1 0 M A « 0 . 

0*1. ~ “ " 

P«P*2. 
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FP«T 
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A4K- A4K*4. _ 
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?0 CfluTIMIP 
C 

WKITF (6,99) (T|«)»«-l,n 

99 HJKmAT ( 1 PH (>15.7/ HD 1 5.7) 

— 12*1 + 1 

H-Fp-T 

IP ( OArtS ( T ) .OT • l « IJ- lo ) H * K / 1 
F*i I w* f)A H S ( P ) 

AnS* T 

90 IP (0AHSU).LF.P.A.s.D.l.Gf.3> r,0 TO 
_j3io. CM^n^up _ 

C " 

1*1-1 

200 InD=16-I 

HkITP ( 6 , 9 A ) P I 'M 
9* FMWnAT ( l P 0 l S . 7 ) 

_ j=* Aiv-i. 

K p TM^tN' 
p»tD 


90S 70? 1 O 


*0670240 
90S /O26o 

w os 7027(7 
*0S702rtO 
*0570290 

90S 703 10 
90-57o32(t 
‘•'OS 7033 m 

90S 703 So 
90S 70 3 60 
ho S 703 Vo 
9 0S 703 Ho 
90S 70 390 
hOS 70400 
H0S704 1 o 


90S? 042n 


90S 704 7m 

200 

90S70490 
90S 70S( If' 
90S TOSlO 


90S 70S3O 
90 s 70S* i r 
*0*>705So 


SIHPTC Pi^C LIST 

C PMi''CT iMivi Suhb^OGKAm AS kPOuIkPU BY I^TPGK AT Ium k uiTImP 

C THIS •rnOTHrf ^PTOK^S COt'»TKML TO T*F I-»TPGkaTIM9 SMHKDUTIwc 

dumhlp pkpcisimv PM^cr i UN pncixj 

DMIlhLP PPPCISIIJ-V X 

cu^" izLUMUi/ )U vi* y t si ( * ) » * »» » xuh f k , r x , i , ph i , j , a x 

2 PKC-x*«(iim*{Ax-M/(i,Himn*n .no-x)**(«x-i . ) 

4 HPTUK-'* 
poll 


SPrV'TKY 


241 

?42 

243 

?44 

24*5 

246 

?67 

24 H 
249^ 

>S0 

?Sl 

2S? 

?63 

2S4 

_2 5S 

25 6 
2S7 
26* 
25 9 
2 60 
2 61 
26? 

263 

264 

265 

266 
2 67_ 
263 ‘ 
269 

2 70 
271 
27? 

2 73 
?74 
2 VS 
276 
2 77 
27 h 
27 9_. 
2 PD ’ 
2*1 
2 P 7 
2*3 
?H4 
?HS_ 
2 H 6 
2 ** 7 

?hH 

2*9 

290 

2.91. . 

29? 

293 

7 9 4 

296 
2 v6 

297 
25* 

299 

300 
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HKdRf <A-« ? M FTm Kk I UK * uFA^ fc-'ST I ■* ATt=‘* 


* 1 h F T C **FAr« K*> 

HAKttMFfFKS IJ^ *FfA KKlUK A, 3 
a ■ -411/ ( l— U ) 
rt ■ KUUf 

I, K .4g4M f {)K AVF^AGf: VAU»h IlF ThK KK' '*dH I L I T Y T HA I A SYSTF-* *UK3> 
ALL UThFK VAKTAHLFS A-<F ►>»-»- I eM -1 1 "> AS I ■ KKHiiKA*. \ 

CdMnO'J / InPuT/U I V ) » V ( b ) t X< { 4 ) ,OPKr K t XviS'FK. I X , I * KH I » J 
h *><")’ FllM m A T U? ) 

VOO FflK-^AI (AFlO.A) 

SOI rMK^AT { S P S ♦ ? ) 

VO? ^iHmaT (4FS.3) 

V03 FHKt* aT { P'S* l > 

H 70 FOKcAT (Fl)*4 f ;MbX,Fb.3»t2<bXtFin.3) f sX,hLi.H} 

KFAO II 

K-AO I S * H SO ) JJ 
KFAO I S ♦ HsU I << 

*FAO i St VOO ) <l>{ I )♦ 1 = 1*11 ) 

kf An (StVni) < v ( j ) ,J*UJJ) 

KPAU (s t vo?) < X <(<)*<= l t <<) 

KF AO (S,V03) rx 

■"d 

I h*0 

1)0 SOU <a\,KK 

on svv Id ♦ II 
on bvH JddJ 
t -STciM I ) *{ I t+X< K) ) 
l»pPhK*Ai-i I iv U 1 • * T U S I ) 

?S XLOWKKaOl I ) V { t ,-X< ( < ) ) 

CALL OKT i 'S IXfl'Mflr) 

Ad u(I )*TX )/ ( I.O-IM I ) ) 

wKITF {4,s70> iMl)fV(J>*X<K)tAtTX f PrtI 

s v* coot I 'OiF 

S^ V C ( T I S'O IF 

son c* ,|M T I '-O *F 
CALL - X I T 
F*d 



3 Ml. 
'A 1 1 > 
3<'3 

3 I * 

3 M s 
304 
3 0/ 
30* 
3rw 

3 in 
311 
3 W 
313 
3 1 4 
31s 
31 s 
31 7 
3\ 3 
3 1 s 
3?.| 

3?1 
33 2 
3 33 
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3?s 
3?7 
3 2 -1 
3? V 
3 3 M 
33 \ 
3 3,> 

3 33 
3 3 + 
33 s 

33 S 
3 3 v 
3^3 
3 -IV 
3 i* i 
3M 
3 u y 
3 A3 
3A/+ 

34 s 
3 A ■■> 
34 ( 
34 3 

34 - 

3-h] 

3s l 
3S? 
3 S3 

35 3 

3bs 

3 s S 
3s 7 
3 s ^ 
3 s V 
3S) 


C 


I ivi f FCKA i' I i Kdliri Jfc -ULL'JW.S n S art*iV»" 
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3*1 

3*2 

3*3 

3*4 

5IKFTC F«'<C LIST 3*5 

f HOC f I r ))M SOKKKIJOKA-'' -e*=MiJ I Kffo KY I v T H GK rt T 1 1 1 xl KMUT I*«l= 3** 

"THIS KflilT I nii= KFTUK'SS COnTkOL TO T> I a'TFCKAT rw SOHKol I \ I of 3*7 

DOOkL* PKFCISIOn Fl^Cfli^ rwC I X ) 3** 

nniJH^H PRECISION X 3*V 

COMMON/ lM>1)T/0( S ) t Vh)fX<|A) ^l^^tXUH^fTX f 1 tHHl , J* ftX 370 

F»\<C« ( ( ( HU ) *AX )-l .0 )/U .0-11 < I J ) ) *{ 1 .00-* ) ** ( AX- I .0 J 371 

4 R F T 1 1 < w Q 7? 

Fun 373 

3 74 
37 3 

F im T * Y "FA N 3 7 * 

377 

37M 

37S 

3^0 

3M 

program 3 prior, -«trf ksti^at-o 3h? 

3*3 

3^4 

IhFTC KOI i T- ? 3*5 

3F* 
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3S0 
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R F A O <5,900) JJ 3s 7 
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M= 1 401 

IF-0 

nri 599 K* 1 « KK 4f ‘3 

H= 1 .0 

no *oo J* l , J J Ui "> 

Call twt (w t H f iw f iF) ^o* 

wkitf (*,9lo> v ( J ) , X< ( < ) f h, ThFT a 407 

*00 Cl i 1 vt j i mi jf 40 3 
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CALL FX IT A l*> 

FiMO 
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ilHHTC AUX 4?l 

SUMKHMTI^ ailX 4?? 
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HKITF (6,80 l ) Y 436 

IF ( Y.FU.O. ) T,n TO 3 437 

THFTA«V( J)-(7/Y) /,- M 

WKITF (6,803) THtTA,Z,Y 

GO TO 4 44 0 

3_^I_JF (6, HO?) 4 M 

6 THFTA - V(\|)-l.O*«l 447 

C ARTIFICIAL VAl.UF I^TcOOtO Til OKIVF KUO T F I Oi 18 K IOHARO A km >1 44 3 

W TO 4 444 

10 THFTA*- 1 , 0 - 70 . 0 #K r 443 

4 FkT-THFTA 4 aS 
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^ ' ,Tr ” <.',i 

44 ^ 
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46 l 
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imYU OAoo (rA8 OOoFi) ROOT I*F ) FOLLOWS 443 
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PURPOSE: To specify the parameters of a prior distribution for two cases: the 

reliability of a unit that either performs satisfactorily throughout a mis- 
sion or does not, and the failure rate of a unit that fails according to the 
exponential distribution. Prediction of demand for spares is considered in 
each case. 

SCOPE : Bayesian analysis is an obvious approach to estimating reliability param- 

eters from such mixed data sources as test results, information on analogous 
components, and engineering estimates. The prior distribution, of necessity 
subjective, is ideally based solely on the latter two sources. This study 
gives a method for specifying the parameters of the prior distribution, re- 
quiring only information corresponding to the most likely value of reliability, 
and to the subjective odds that the error in this estimate is less than a given 
percentage. Testing data are then merged with the prior distribution via Bayes 
rule to obtain a posterior distribution. Roughly speaking, the spread of the 
prior distribution is inversely proportional to the degree of prior belief and 
determines how heavily the distribution will be weighted when it is combined 
with test data. Tables of the parameters of the prior distribution, computed 
according to the method described, appear in the Appendix. 

BACKGROUND : This research was done by RAND for the National Aeronautics and Space 

Administration in connection with the Apollo mission reliability assessment 
study. It should be of interest to those working on reliability estimation, 
allocation of investment among system components to achieve maximum system 
reliability, and stockage applications. 





